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Co0 2 -Layer-Thickness Dependence of Magnetic Properties and Possible Two 
Different Superconducting States in Na^CoC^ • 2/H 2 
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In order to understand the experimentally proposed phase diagrams of Na^Co02 • 3/H2O, we 
theoretically study the CoCh-layer-thickness dependence of magnetic and superconducting (SC) 
properties by analyzing a multiorbital Hubbard model using the random phase approximation. 
When the Co valence s is +3.4, we show that the magnetic fluctuation exhibits strong layer- 
thickness dependence where it is enhanced at finite (zero) momentum in the thicker (thinner) 
layer system. A magnetic order phase appears sandwiched by two SC phases, consistent with 
the experiments. These two SC phases have different pairing states where one is the singlet 
extended s-wave state and the other is the triplet p-wave state. On the other hand, only a 
triplet p-wave SC phase with dome-shaped behavior of T c is predicted when s=+3.5, which 
is also consistent with the experiments. Controversial experimental results on the magnetic 
properties are also discussed. 
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A possible unconventional superconductivity (SC), 
particularly a spin-fluctuation-mediated one, has been 
expected in Na^CoO^ • yH^O. 1,2 Hence, its normal-state 
magnetic property has been intensively studied, although 
the results are rather scattered. 2-10 

The bulk susceptibility x shows Curie- Weiss type up- 
turn below ^130 K with decreasing temperature. Several 
groups have observed a similar upturn in Knight shift 
(K) and a linear K-\ plot in NMR and /xSR measure- 
ments, suggesting an increase of spin fluctuation at/near 
q = Q 2,5,6,9 Among them, Ishida et al. studied the re- 
lation between \/T\T and X- 3,6 They first reported an 
identical T dependence of 1/T\T and \ up to T c , which 
suggests a dominant ferromagnetic (FM) fluctuation at 
q~0. 3 Later, they reported a slightly different behavior 
in another sample where the upturn of x 1S weaker than 
that of 1/7]T, which suggests dominant spin fluctuations 
at q ~ but not at q = 0. 6 On the other hand, Ning et 
al. and Mukhamedshin et al. observed T-independent be- 
havior of K despite the strong T dependence of \jT\T, 
which suggests dominant spin fluctuations at g ^ 0. 7 ' 8 
We also note that a neutron-scattering experiment did 
not detect any evidence for spin fluctuations. 10 

On the other hand, a relationship between T c and 
Co02-layer thickness has been pointed out by several 
groups. 11-18 In particular, Sakurai et al. determined x- 
T phase diagrams where x is the Na content. 11 ' 12 Here, 
x scales with the Co02-layer thickness, i.e., a larger- x 
sample has thicker C0O2 layers. At the same time, they 
found that the Co valence, s, is constant at ^+3.4 al- 
though x changes because of the presence of H30 + ions 
in Na-layers. As a function of x, they found successive 
three phases of SC (SCI), a magnetic order (MO) and 
another SC (SC2). Note that s is directly related to the 
number of i2g-electrons per Co ion (ji t 2 g ) as n t 2 g = 9 — s. 
On the other hand, for samples with a slightly diffcr- 

*E-mail address: mochizuki@rikon.jp 



ent value of s (s ~+3.5), only one SC phase appears, 
and T c shows a dome-shaped behavior as a function of 
x. This indicates that a subtle change in the lattice pa- 
rameter and that in the Co valence affect drastically the 
electronic properties. 

Motivated by these findings, we previously studied ef- 
fects of CoOg distortion on the band structures. 19 We 
constructed an eleven-band tight-binding (TB) model 
including the Co 3d and the O 2p orbitals, which re- 
produces very well the LDA data for the bilayer-hydrate 
system of ref. 20. In the case of s=+3.4, we found that; 
(i) Fermi surface (FS) with double ai g -band cylinders 
around the T point is realized in a system with thick 
C0O2 layers. This FS was referred to as FS1 (see the bot- 
tom figure of Fig. 1(b)). (ii) FS with a single ai g -band 
cylinder and six e' ff -band hole pockets (FS2) is realized in 
a thin-layer system, (hi) In the moderate-thickness case, 
another type of FS with double a\ g cylinders and six 
e' g pockets (FS3) is expected. We discussed that this FS 
deformation can explain the experimental s=+3.4 phase 
diagram with three successive phases. 19 

In this letter, we perform microscopic calculations on 
the SC gap structures as well as the nature of spin 
fluctuation since their knowledge is essentially impor- 
tant to understand the experimental results. We resolve 
the above discrepancies of the experimental results of 
the character of magnetic fluctuation. The Co02-layer- 
thickness dependence of magnetic fluctuation and SC 
states are studied by constructing the multiorbital Hub- 
bard model with threefold Co t2 g orbitals and by ap- 
plying the random phase approximation (RPA). When 
s is +3.4, we show that the spin fluctuation is critically 
enhanced with decreasing T toward MO in the moderate- 
thickness system with FS3 (see Fig. 1(b)). This is in 
agreement with the existence of MO phase in the exper- 
imental phase diagram. By solving the Eliashberg equa- 
tion, we show that the singlet extended s-wave pairing is 
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expected in the thick-layer systems with FS1, while the 
triplet p-wave pairing is expected in the thin-layer sys- 
tems with FS2. For the s=+3.5 case (see Fig. 1(c)), in 
contrast, the magnetic instability hardly occurs and only 
a triplet p-wave SC phase with dome-shaped T c behavior 
appears, which is also consistent with the experimental 
s=+3.5 phase diagram. 

First, let us discuss the TB model used in this paper. 
As noted before, we developed the eleven-band TB model 
in ref. 19. But that model is not covenient for numerical 
calculations because of the large degrees of freedom. In- 
stead, we construct a simpler three-band TB model with 
only Co tig orbitals, which reproduces the layer-thickness 
dependence predicted in the previous eleven-band anal- 
ysis. 

The obtained three-band TB Hamiltonian is given by 

#3TB = Ek,m,n,, £ k ^L/kna with 

e 77 = 2ii cos fc 77 + 2* 2 cos /c 77 + 2* 2 cos^ 77 + W 

+ 2* 5 cos^ 77 + A: 77 ) + 2* 5 cos(fc 77 - fc 77 ) 

+ 2* 6 cos(/fc 77 + 2fc 77 ) + 2* 9 cos(2fc 77 ) 

+ 2< 10 cos(2fc 77 ) + 2* 10 cos(2fc 77 + 2fc 77 ), 

eg = 2*3 cos kg + 2*4 cos(/c 77 ' + fc 77 ' ) + 2* 4 cos kji) 

+ 2< 7 cos(/c 77 ' + 2kg ) + 2* 7 cos(fc 77 ' - fc 77 ' ) 

+ 2< 8 cos(2fc 77 ' + kg) + 2t n cos(2fc 77 ') 

+ 2*i 2 cos(2fc 77 ' + 2kg) + 2< 12 cos(2fc 77 ') - A/3. 

Here, 7 and 7' represent xy, yz and zx orbitals, and A 
denotes the trigonal crystal field (CF) from the O ions. 

For definitions of kg ' and /c 77< , one could see ref. 25. 

Compared with the previous TB model in ref. 25, the 
present TB model is much more elaborate including ad- 
ditional transfer integrals. It reproduces subtle features 
of LDA results of ref. 20. Note that the numbering of 
transfer integrals slightly differs from ref. 25. 

In Fig. 1, we show (a) band dispersions and (b) FSs for 
several values of OCoO-bond angles (<^ocoo)- Here, the 
angle 0oCoO expresses the Co0 2 -layer thickness where 
a thinner C0O2 layer has a larger 0ocoO value (see 
Fig.l (a) in Ref. 19). The LDA data used in the pre- 
vious study was calculated using experimental structure 
data with </>oCoO = 97.5°. The parameter values of the 
present model are (*i, t 2 , *i 2 , A) = (35.0, —22.0, 153.5, 
46.1, -17.7, -14.9, 3.10, -52.4, -41.0, -27.6, 8.16, 4.98, 
80.0) for the case of 0oCoO = 97.5° where the unit is 
meV. Note that FS1 (FS2) is reproduced for the thick 
(thin) layer system, while FS3 is reproduced for the mod- 
erate case. This type of FS-topology variation is referred 
to as Case C in the previous study. 19 

By further adding the Coulomb interaction term Hi nt ., 
we obtain the multiorbital Hubbard model; H mo = 
i?3TB+#int. where H int . = Hu + Hu'+Hj^+Hj> + H v . 
As in the previous studies, 25,26 the terms Hjj and Hw 
represent the intra- and inter-orbital Coulomb inter- 
actions, respectively, and Hj H and Hj> represent the 
Hund's-rule coupling and the pair hopping, respectively. 
These interactions are expressed using Kanamori param- 




Fig. 1. (a) Band dispersions for various 0oCoO values calcu- 
lated from the three- band tight-binding model -H3TB ■ Horizontal 
dashed lines denote the Fermi level for s=+3.4. (b) Deformation 
of the Fermi surface with varying </>oCoO f° r the s=+3.4 case, 
(c) Deformation of the Fermi surface for the s=+3.5 case. 



eters, U, U' , Jh and J', which satisfy the relations; 
U 1 = U—2Jn and Jh = J'- In this paper, we include the 
last term Hy — V ^ UiTij representing the Coulomb 
repulsion between adjacent i and j sites. 

We analyze this model by applying RPA. In the 
present three-orbital case, the Green's function G is ex- 
pressed in the 3x3-matrix form corresponding to the xy, 
yz and zx orbitals. The irreducible susceptibility \° has 
a 9x9-matrix form. The singlet (triplet) pairing inter- 
action r s (r*) is expressed using the interaction matri- 
ces. For detailed expressions of G, x°, T s and T*, see 
ref. 25. Note that the matrices U s (q) and U c (q) slightly 
differ from those in ref. 25 since the present model in- 
cludes long-range Coulomb repulsion Hy- The matrix 
elements U^Jq) (U^Jq)) are U (U + 2V(q)) for 
m = n = = v, J H (2U' + 2V(q) - Jh) for m = 
n 7^ [i = v, U' {—V + 2Jh) for m = (a ^ n = is, 
J' (J') for m = v 7^ n = fj,, and for others, where 
V(q) = 2V r [cos(qi) + cos(q 2 ) + cos(gi + q 2 )] with q x = 
\/3/2q x — l/2q y and q 2 — q y . We discuss the nature of 
SC by solving the Eliashberg equation. Calculations are 
numerically carried out with 128x128 fc-meshes in the 
first Brillouin zone, and 1024 Matsubara frequencies. 

We first discuss the results for the s=+3.4 case, which 
are calculated taking [7=0.5 eV and Jh=0.05 eV. From 
now on, the energy unit is eV. In Fig. 2(a), we display the 
spin susceptibility Xs(<?) for several values of T and the 
angle </>oCoOj which shows strong layer-thickness depen- 
dence. For the thick-layer cases with FS1 (</>oCoO = 96.5° 
and 97.0°), Xs(l) nas peak structures at finite momentum 
q = Qs in the low-T region. This q = Q s fluctuation is 
induced by the electron scattering between the inner and 
outer dig FSs owing to the intra-orbital Coulomb repul- 
sion U. Actually, Q s is the wave number which bridges 
these inner and outer FSs. On the other hand, for the 
thin-layer cases with FS2 (0ocoO = 98.0° and 98.5°), 
Xs(q) has a FM peak at q = 0. This FM fluctuation 
is induced by the inter-orbital Hund's-rule coupling Jh 
between a\ g and e' g FSs as discussed previously. 25 
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Calculated results on the spin susceptibility Xs(l) for the 
s=+3.4 case, (a) Momentum dependence of Xs(l) for several T 
and </>oCoO values, (b) Xs(?) a t 9=0 (Xs(? = 0)) and maximum 
Xs(q) (max[x s (g)]) plotted as functions of T for several 0oCoO 
values. 



For the moderate layer-thickness systems with FS3 
(<^OCoO = 97.0°-98.0°), a critical enhancement of Xs{l) 
with decreasing T occurs, which is consistent with the ex- 
istence of MO phase sandwiched by two SC phases. 11,12 
Indeed, Xs(l) for 0ocoO = 98.0° diverges at k^T = 
0.012. This magnetic instability is caused by cooperative 
contributions from the intra-band scattering between the 
inner and outer a\ g FSs and the inter-band scattering be- 
tween the aig and e' g FSs owing to the FS3 geometry. In 
addition, the structure of density of states (DOS), where 
both aig and e' g orbital components are large, is also re- 
sponsible for the magnetic instability. 19 Indeed, as will 
be shown later, magnetic instability does not occur for 
the s=+3.5 case with FS3 because of a small a± g orbital 
component of DOS. 

The puzzle in the NMR/NQR and fiSR results can 
be solved by considering layer-thickness dependence of 
Xs{q)- The quantities K and 1/T\T scale, respectively, 
with Xs(<Z = 0) and maximum Xs(<l) ( max [Xs(9)])- These 
quantities obtained in the present theory are plotted 
in Fig. 2(b). In the thick-layer case (</>oCoO = 96.5°) 
with FS1, Xs(q — 0) is rather T-independent since the 
spin fluctuation is not FM, which is consistent with 
the reports about T-independent NMR K from Ning et 
al and Mukhamedshin et al. 7 ' s In the thin-layer case 
(<^OCoO = 98.5°) with FS2, on the other hand, the spin 
fluctuation is FM, and Xs(l = 0) and max[x s (<z)] show 
an identical increase, which reproduces the x-(l/TiT) 
relation in ref. 3. Finally, in the moderately-thick case 
(^OCoO = 97.0°), Xs(q = 0) shows weaker increase than 
max[x s (g)], which reproduces the x~{^/T\T) relation in 
ref. 6. Indeed, the sample in ref. 6 turned out to have 
moderately thick layers according to the measured NQR 
frequency vq, which is consistent with the present result. 

Next, we discuss the SC properties. In Fig. 3, T c 
for several pairing states are plotted as functions of 
V both for the thick-layer case with 0oCoO = 96.5° 
(Fig. 3(a)) and for the thin-layer case with </>ocoO — 
98.5° (Fig. 3(b)). We find that in the thick-layer system 
with FS1, the singlet extended s-wave state is stabilized 
in the wide range of V value while the triplet p-wave 




0.15 



extendeds-wave pr-wave 



py-wave 



(f) 



1.4 



0.6 



0.2 



•s=+3.4(«t2 g =5.6) 
<j)OCoO=96.5° 
teJ=0.012 eV 
F=0.15 eV 

0.10 eV 
0.05 eV 

0.0 eV 




K \ 





(g) 



5=+3.4(«t2g=5.6) 

Magnetic i 




K 



M K 



thicker < C0O2 layer 

Na content* 



triplet"S€ 
(p-wave) 



<J>OCoO 
-> thinner 



Fig. 3. Calculated results on superconducting properties for the 
s=+3.4 case, (a) T c for several pairing states plotted against V 
for a thick-layer case with </>oCoO = 96.5°. (b) that for a thin- 
layer case with </>oCoO = 98.5°. (c) Gap structure of the extended 
s-wave pairing on the FS1 for </>oCoO = 96.5°. (d) and (e) Gap 
structures of the p x and p y pairings on the FS2 for </>oCoO = 
98.5°. (f) Charge susceptibility Xc{q) in the momentum space 
for several V values, (g) Schematic T-c^oCoO phase diagram for 
the s=+3.4 case. 



state is stabilized in the thin-layer system with FS2. Fig- 
ures 3(c)-(e) show the fc-dependence of the obtained SC 
gaps with (c) extended s-wave, (d) p^-wave, and (e) p y - 
wave symmetries. This result shows that two SC states 
with different symmetries arc possibly realized in this 
material depending on the Co02-layer thickness. 

The extended s-wave gap here obtained for the FS1 
case is equivalent to the one previously proposed by 
Kuroki et al. 21 The signs in this gap are the same within 
each FS but are opposite between the inner and outer 
FSs. This gap structure is stabilized not only by the 
spin fluctuation \ s (?) but also by the charge fluctuation 
X°((?) as discussed previously. In fact, when V = 0, A 
for p-wave state is slightly larger. The extended s-wave 
state dominates when we introduce small V (at least 
V/U ~ 0.02). Note that in the expression of the singlet- 
pairing interaction, the contribution from X s (q) an d that 
from x c (l) have different signs. The repulsive contribu- 
tion from X s (?) around q — Q s favors the sign change 
between the inner and outer FSs. On the other hand, 
the attractive one from x c (l) favors the same sign in the 
outer FS. This is because x c (g) tends to have peak struc- 
tures around the M-points as well as near the K-points 
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Fig. 4. Calculated results for the S=+3.5 case, (a) Inverse of 
maximum \k (</) plotted against T for several flioCoO values, (b) 
T-ipoCoO phase diagram, which shows a dome-shaped T c of the 
triplet p-wave pairing state. 



when V is introduced (see Fig. 3(f)). These wave num- 
bers correspond to those across the outer FS. 

As for the p-wave state, which is realized on FS2, the 
p x and p y states are degenerate in the triangular lattice. 
Below T c , a linear combination of these two should be 
realized. 25,26 As shown in Figs. 3(d) and (e), the gap am- 
plitude is large on the a\ g band while it is markedly small 
on the e' g band. This is in contrast to the previous results, 
which show dominant gaps on the e' g pockets. 25-28 This 
difference is caused by the difference of TB models used. 
In the previous model, the a\ g band has a steep slope 
at the Fermi level resulting in a small DOS, which is 
unfavorable for the gap opening on the a\ g FS. On the 
other hand, the present model has the a\ g band whose 
slope is rather gradual, resulting in a larger DOS and the 
dominant SC gap on the a\ g FS. Since the present model 
reproduces the LDA band structure much more precisely, 
we consider that the present result is more realistic. 

Some groups observed decreasing K below T c in the 
NMR experiments. 22-24 We speculate that samples used 
by Kobayashi et al. are located in the predicted singlet 
extended s-wave phase according to the measured NQR 
frequency vq. 29 On the other hand, T- independent be- 
havior of K below T c may be observed if we measure a 
sample with the triplet p-wave phase. However, synthesis 
of such a sample is rather difficult and should be done 
carefully since the triplet SC is quite fragile against im- 
purities and oxygen defects. There exist discrepancies in 
the experimental data on H C 2 and specific heat. 2 These 
can be also resolved if we consider the two different SC 
states with different FS topologies. The inner a\ g FS in 
FS1 and the e' g pockets in FS2, which reproduce the ex- 
perimental data, are proved to be quite small. These is- 
sues will be discussed in detail elsewhere. The predicted 
FS deformation should be detected in the bilayer-hydrate 
materials by the bulk-sensitive angle resolved photocmis- 
sion spectroscopy. 30 

Finally, let us discuss briefly the case with s=+3.5. In 
this case, the magnetic instability is weak as compared 
to the s=+3.4 case. In Fig. 4 (a), we show l/max[x s (g)] 
plotted against T for several </>oCoO values, which are cal- 
culated taking [7=0.58, J H =0.05 and V=0.0. This figure 
shows no critical enhancement of Xs(q) even though we 
use slightly larger Coulomb parameters than those used 
in the calculations for the s=+3.4 case. As for the SC 



state, it is proved that only the triplet p-wave pairing 
is stabilized. As shown in Fig. 4 (b), T c shows a dome- 
shaped behavior as a function of 0ocoO > which is consis- 
tent with the experiments. 11, 12 

To summarize, we have studied the Co02-layer- 
thickness dependence of magnetic and SC properties in 
NazCo02 • yEbO. By analyzing the multiorbital Hub- 
bard model using RPA, we have reproduced the ex- 
perimentally obtained s =+3.4 phase digram contain- 
ing successive SCI, MO and SC2 phases as well as the 
s =+3.5 phase diagram containing one SC phase with 
dome-shaped T c behavior. We have shown that two SC 
phases for s =+3.4 have different pairing states where 
one is the singlet extended s-wave state and another is 
the triplet p-wave state, while the SC phase for s =+3.5 
has the p-wave state. We also discuss that the puzzling 
NMR/NQR and pSR results on the character of mag- 
netic fluctuation can be understood by considering the 
strong layer-thickness dependence of the magnetic fluc- 
tuation. 31 
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